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We use linear response theory to study the transverse force generated by an external electric field 
and hence possible charge Hall effect in spin-orbit coupled systems. In addition to the Lorentz force 
that is parallel to the electric field, we find that the transverse force perpendicular to the applied 
electric field may not vanish in a system with an anisotropic energy dispersion. Surprisingly, in 
contrast to the previous results, the transverse force generated by the electric field does not depend 
on the spin current, but in general, it is related to the second derivative of energy dispersion only. 
The transverse force always vanishes in the system with an isotropic energy dispersion. However, 
the transverse force may also vanish in some systems with an anisotropic energy dispersion such as 
the 2D k-cubic Dresselhaus system. Furthermore, we find that the transverse force does not vanish 
in the Rashba-Dresselhaus system. Therefore, the non-vanishing transverse force acts as a driving 
force and results in charge imbalance at the edges of the sample. This implies that a non-zero Hall 
voltage can be detected in the absence of an external magnetic field in anisotropic systems such as 
the Rashba-Dresselhaus system. The estimated ratio of the Hall voltage to the longitudinal voltage 
is ~ 10"'^. The disorder effect is also considered in the study of the Rashba-Dresselhaus system. We 
find that the transverse force vanishes in the presence of impurities in this system because the vertex 
correction and the anomalous velocity of the electron accidently cancel each other. Nonetheless, we 
believe that the transverse charge imbalance can be detected in the ballistic region by measuring 
the Hall voltage. Our interesting prediction would stimulate measurements of the Hall voltage in 
such spin-orbit coupled systems with an anisotropic dispersion as the Rashba-Dresselhaus system 
in the near future. 

PACS numbers: 71.70.Ej, 71.55.Jv, 72.10.-d, 73. 23. Ad 



I. INTRODUCTION 

The equation of motion of a charged particle in the 
presence of an electromagnetic field is determined by the 
classical Lorentz force |H. When only an electric field is 
present, the charged particle is accelerated along the elec- 
tric field (longitudinal motion). In the presence of only 
a magnetic field, the charged particle is deflected by the 
magnetic field when the particle is in motion (cyclotron 
motion) . In the 2D plane, the combination of longitudi- 
nal and cyclotron motion would result in the appearance 
of the classical charge-Hall efl^ect [2] when the applied 
electric field and magnetic field are parallel and perpen- 
dicular to the plane, respectively. In the absence of a 
magnetic field, the charge Hall effect would disappear in 
the classical regime. 

When the spin degree of freedom is considered, the 
charged-particle motion further changes; in particular, 
the orbital motion now depends on the spin direction 
(spin-orbit interaction). The anomalous motion of par- 
ticles that is induced by spin-orbit interaction has at- 
tracted much attention. The semiclassical force exerted 
on a spin- 1/2 electron in the presence of an electromag- 



* Electronic address: 
t Electronic address: 



twchen@phys.ntu.edu.tw| 



gyguo®phys.ntu.edu.tw| 



netic field has been investigated Q. The force acts on 
the spin according to SU(2) non-abelian gauge theory 
because the spin-orbit interaction (cr x E) • p, which is 
taken to be the kinetic momentum, plays the role of the 
gauge field It has been shown that the equation of 
motion of an charged particle in a crystal environment 
can be altered by the Berry curvature of band struc- 
ture 0, Q when the spin-orbit interaction is considered: 

X = ^l^g^^ — k X B„ (k) . The Berry curvature correction 
(— kx B„(k)) can also be derived by using the Feynmann 
path integral method . The Berry curvature plays the 
role of an effective magnetic field and thus results in the 
anomalous deviation of the classical trajectory. In this 
context, an interesting phenomenon was observed in the 
study described in Ref. in the adiabatic approxima- 
tion, the transverse velocity due to Berry curvature ex- 
ists even in the absence of an external magnetic field. It 
was further shown that in the adiabatic sense, the non- 
vanishing Berry curvature near the Brillouin zone center 
of a hole-doped (p-type) semiconductor would lead to in- 
trinsic spin-Hall effect Q . This discovery led to an inter- 
esting result: the magnetic monopole would exist in the 
crystal momentum space of solids Furthermore, the 
concept of the anomalous velocity due to Berry curvature 
has been used in ab initio relativistic band structure cal- 
culations of the intrinsic spin Hall conductivity in p-type 
semiconductors 9] and also Ft metal [lo| . 
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Due to the existence of Berry curvature , the intrin- 
sic force of the k-linear Rashba system is studied in semi- 
conductor wires in ballistic regime [lij. The force was 
semiclassically defined as the time derivative of the ki- 
netic momentum in the Heisenberg picture. It was found 
that the deflection of the center of spin-polarized wave 
packet in the transverse direction can be explained by the 
existence of the intrinsic force [ll[ . The concept of intrin- 
sic force due to the spin-orbit interaction in semiconduc- 
tor is further generalized to the non-relativistic limit of 
the Dirac equation in which the spin-orbit interaction 
is still linear in momentum. In the non-relativistic limit 
of the Dirac equation, it has been shown that a semiclassi- 
cal force is exerted on the moving particle with spin ■ 
For the 2D system, the effective Hamiltonian describ- 
ing the Rashba-Dresselhaus system is found to have the 
structure of the SU(2)x U(l) gauge field The trans- 
verse force derived from the corresponding four-force in 
the Rashba-Dresselhaus system agrees with the result in 
Ref. • Recently, the inverse spin-Hall effect driven by 
the spin motive force in a ferromagnetic conductor has 
been studied [ij . The Hall angle in this case is shown to 
be greater than that in the anomalous Hall effect since 
skew and side jump scattering are considered pl| . 

However, the force in the linear response to an exter- 
nal electric field in a spin-orbit coupled system has not 
yet been studied. This is very important in the study 
of forces generated by external fields because the force 
operator is significant only when the expectation value 
is taken into account. In the present paper, we study 
the force in spin-orbit coupled systems with the generic 
two-band effective Hamiltonian [Eq. ([T])] by using linear 
response theory. The force is semiclassically defined as 
the time derivative of the kinetic momentum mv. The 
velocity v can be separated into two terms, v = + v^, 
where = p/m and is the anomalous velocity due 
to the spin-orbit interaction. The spin-orbit interaction 
induced force (SOIF) originates from the time deriva- 
tive of mVa- In this study, we shall focus on the SOIF 
in the linear response to an external electric field. Our 
main result is that the SOIF depends only on the second 
derivative of energy dispersion and not on the spin Q 
or spin current 12]. We provide an analytic expression 
for the SOIF and show that the non- vanishing transverse 
SOIF may exist when the system is not rotationally in- 
variant. This also implies that the charge imbalance at 
the edges of the sample may occur in a spin-orbit cou- 
pled system because the non-vanishing transverse force 
acts as a driving force in the system. 

Our present paper is organized as follows. In Sec. II, 
we derive the force from the generic effective Hamiltonian 
by using the semiclassical definition, that is, the time 
derivative of kinetic momentum mv. In Sec. Ill, we use 
linear response theory to investigate the properties of the 
force under a weak and constant electric field. In Sec. 
IV, we study the Rashba-Dresselhaus system and also 
consider the presence of weak disorder. Our conclusions 
are presented in Sec. V. 



II. FORCE GENERATED BY SPIN-ORBIT 
INTERACTION 

In this section, we will derive the intrinsic force in the 
linear response to an electric field. We assume that the 
effective Hamiltonian can be partitioned into two terms: 
a kinetic energy term and spin-orbit interaction term. 
In the absence of spin-orbit interaction, the Hamiltonian 
would simplify to the free electron (hole) model. In this 
case, we consider the following effective Hamiltonian that 
includes the spin-orbit interaction: 

Hq = —— + hso, (1) 
Zra 

where the first term of Eq. ^ represents the kinetic en- 
ergy for a free effective mass m and hgo represents the 
spin-orbit interaction. We do not consider the external 
magnetic field, and thus, the time reversal symmetry is 
preserved. By spectral decomposition, the spin-orbit in- 
teraction hso can be written as 

hso = A„(k)7'„fc 

(2) 

Vnk = \nk){nk\, 

against the background that the eigenenergy of the 
Hamiltonian is i?„k — ^^-^ + A„(k); in other words, 
Ho\n\i) — En]i\nk). Vnk is the projection operator 
that satisfies the relation Vnk'Pmk = VnkSmn- The in- 
dex n denotes the band index, and A„(k) is the en- 
ergy dispersion of spin-splitting for the n-th band. It 
must be emphasized that the dispersion A„(k) is not 
limited to the odd function of k in which the system 
lacks an inversion center. For GaAs, the parity selec- 
tion rule does not allow the appearance of a fc linear 
term for the hole band, and the spin-orbit interaction 
term of the Luttinger Hamiltonian is quadratic in k. For 
a spherical approximation of the Luttinger Hamiltonian 
(see Ref. [7]) Ho = + ll2)k^ - 272 (k • S)^], we 

have hso — 2^'^2[|fc^ — 2(k • S)^] for the spin-orbit in- 
teraction, m = mo/71 for the free effective mass, and 
A±3(k) = -272!^ and A±i(k) = +272^ represent 
the spin-splittings of heavy and light holes, respectively. 
On the other hand, for a 2D system, the spin splitting 
can be simply written as A„(k) = — nA(k), which indi- 
cates that the magnitude of spin-splitting in each band 
is the same. In the following discussions, we will show 
that only a fc-linear term can reproduce the result given 
in Ref. [12l |. For terms containing higher orders of fc, in 
general, the force operator is not related to the spin cur- 
rent. We will show that in the 2D system, the force is 
actually related by the Berry vector potential. The com- 
mon 2D effective Hamiltonian that includes spin-orbit 
interaction can also be written in the form of Eq. ([1]) 
(see also Ref. [IBl)- It should be pointed out that the 
knowledge of explicit form of the energy dispersion A(k) 
and eigenfunction |nk) is not necessary in the following 



3 



derivations. We will derive an analytic expression for the 
linear response of the semiclassical force. 

In the presence of a constant electric field, the potential 
can be written as — gE • x, and we have 



(3) 



where q is the charge of the particle (g = — |e| for an 
electron). The notation x refers to the gauge invariant 
position operator [l^. The diagonal terms of the gauge 
invariant position operator vanish, and the off-diagonal 
terms is equivalent to the replacement of x with i^. 
The force in the following derivation is defined as the 
time derivative of kinetic momentum mv, and we have 



(4) 



where the velocity operator is defined as v = jf^ [^; H] = 
where the gauge invariant position operator does 
not change the original commutator [l6j . The velocity 
operator is composed of two terms: the free electron 
(hole) velocity v/ = p/m and the anomalous velocity 
Va = dhso/dp^ which is induced by the spin-orbit inter- 
action. By the use of equation v = v/ + v^, the semi- 



classical force F = 
terms: 



dt 



||[v, i/] is composed of two 



F^F^+F^^-^, (5) 

where F^ = j^[vf,H] can be written as (the Einstein 
summation convention is used) 



qSijEj, 



(6) 



which is the Lorentz force. F'^^' = Ifva, i?] is the SOIF 
operator. For convenience in the following discussions, 
the SOIF operator is rewritten as the sum of two terms: 



+ F0. The force F^ = ^[j^,-qF ■ x] de- 
pends on the anisotropic properties of energy dispersion, 
and it can be shown that 



qra d^hgo 
h? dkidkj 



(7) 



The force F*^ originates from the commutator of the 
anomalous component of the velocity and the spin or- 
bit interaction of the Hamiltonian, namely, [vq , hso] , and 
it can be written as 



F,^ 



m dhso 



(8) 



In general, Eq. ([8]) has a complicated form, but in a 2D 
spin-orbit coupled system, it can be shown that the force 
Ff^ [Eq. ([5])] is related to the Berry vector potential. This 
can be seen as follows. For 2D systems, the Hamiltonian 
can generally be written as i/p = %^+^(k)cra: — S(k)crj,, 
with energy dispersion Enk — — + B'^. It can 

be shown that 



rn_ dhso_ , , _ 



86 

dki 



(9) 



and 17 



2A 

h ■ 



where 0(k) is defined as 6* = tan ^ (s^^ 
For an appropriate choice of phase of eigenstate, one can 
obtain (nk|(— i)^|nk 



1 ae 

2 ak- 



in Ref. [17| , it was shown 



that the spin-Hall conductivity can be expressed in terms 
of anti-commutator of Eq. ([5]) and spin current operator. 

On the other hand, it is interesting to note that for 
k-linear systems such as the Rashba-Dresselhaus system, 
Eq. (El) can be further written as ^{a^ - 13'^){3'''^ x e^), 
where J''^ = ^{v, Sx} is the conventional definition of 
spin current, as shown in Ref. p^ . where the spin 
force corresponding to F*^ is derived by considering the 
non-relativistic limit of the Dirac equation. However, 
the force F*^ is not always related to the conventional 
spin current. For example, for the 2D k-cubic Dressel- 
haus system [l5|, [l^ , where the spin-orbit interaction is 
hso = f^Dikxkydx — kyk'^Uy), and the energy dispersion 

is A(k) = Pokkxky, we have dO/dkx = Therefore, 
Eq. ^ can be written as k1ky{kx(Tz), where the 

spin-Hall current can not be determined. 



III. LINEAR RESPONSE TO WEAK ELECTRIC 
FIELD 

Using Ehrenfest's theorem, the time derivative of the 
expectation value of the velocity is equal to the commu- 
tator of velocity and system Hamiltonian, i.e.. 



d_ 

m 



"*^(t)v*(i)] = Re 



*t(t)l[v,i/]vl/(i) 

in 



(10) 



where Re[- ■ ■] represents the real part of [• • • ]. The state 
vector ^{t) satisfies the Shrodinger equation = ih,^ 
and if = iJo - gE • X [Hq is given by Eq. IJ)). In the 
derivation of Eq. pO)l . the fact that {vi,H} is a hermi- 
tian operator is used. Accordingly, the left-hand side is 
defined as the force in the system, and the right-hand 
side of Eq. (fTU)) can be expanded up to the first order of 
the applied electric field. In linear response theory, the 
expectation value of the operator O can be written as 



(^) = yE/»k(nk|0|nk) 



;^^/„k2i?e[(nk|0|,ik)i]+o(i?2), 



(11) 



nk 



where /nk is the Fermi-Dirac distribution and V is the 
volume of the system. In the right-hand side of the equal- 
ity, the first term of Eq. pTjl gives the expectation value 
of the unperturbed wave function and the second term 
of Eq. pT|) represents the linear response to the weak 
electric field. The perturbed wave function to first order 
of perturbation V , denoted as |fik)^, is given by 



ik)i = 



'klylnk) 



En 



En'V 



(12) 
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In our case, the perturbation is V = —qE ■ x. It must be 
emphasized that the gauge invariant position operator 
does not change the matrix element of (n'k|y|nk) be- 
cause only the interband transition is considered in Eq. 

(na) [13. 

In the following calculation, the expectation value is 
evaluated to the first order in the electric field. For the 
Lorentz force [Eq. ©J , we have (to the first order in 
the electric field) 



(13) 



where ric is the carrier concentration. The Lorentz force 
acts only in longitudinal direction. It is important to note 
that the Lorentz force is derived from the commutator of 
the kinetic momentum ^ and potential V ~ — qE • x. In 
single-layer graphene and bilayer graphene, there is no ki- 
netic energy in effective tight-binding Hamiltonians p^ . 
and thus, the charged particle does not experience the 
Lorentz force in either of the systems. In calculating the 
linear response of [Eq. ([7])] and F*^ [Eq. ©J, wc note 
that Eq. ([7]) includes the effect of electric field but Eq. 
([5]) does not. By using Eq. ^ and the first term of Eq. 
PT|) . it can be shown that 



(nk 



qmEj d^hso 
h"^ dkodki 



9^A„(k) qm 



dkjdkj 



Ink) 



where 



2A„i?e 



2Re 



,a(nk) 



id{nk) d{nk) , 
' dh dkj ' 
d{nk) ^ 



dkj 



(15) 



The calculation of second term of Eq. (fTTj) for is not 
necessary because it gives a second order of an electric 
field. The linear response of Eq. ([5]) is obtained from the 
expectation value 2i?e(nk|i[^^, /iso]|nk)^ because the 
first term of Eq. (|lip vanishes, i.e., it can be shown that 
(nk| i [^^, /iso] jnk) = 0. For interband transition, the 
gauge invariant position operator x can be replaced by 
and we have 

2i?e(nk|i[?J^,/i,o]|nk)i 



i dki 

2qEj ^ (A„-Af)i?e 

^2 ^ij^J- 



'^^"^)|^k)(^k|-^^"^)^ 



dki 



dki 



(16) 



We note that the restriction £ n in the first equality of 
Eq. (|16p is not necessary because the term corresponding 
to £ = n is zero. Furthermore, the second equality is 



carried out by using the fact that /iso = 
and Y.I Kk)(^k| = 1. We obtain 



nk 



(17) 



It is interesting to note that Eq. (|17p exactly cancels 
the second term of Eq. This implies that in fact, 

the force (F*^) does not exist in the linear response to a 
static electric field. The resulting mean SOIF (Ff^^) = 
(Ff+Ff) is (the Einstein summation convention is used) 



SOI\ 



{Ff + F[')^qK,,E, 



where 



9A„(k) 
' dkidkj 



(18) 



(19) 



Therefore, in the system with spin-orbit interaction, the 
particle would experience two kinds of forces when only 
an electric field is present: 



(Fi) = q{ncdij + Ktj)Ej, 



(20) 



(14) where Kij is given in Eq. (fTO)) . The first term repre- 
sents the Lorentz force and the second term represents 
the force generated by the spin-orbit interaction. Physi- 
cally, Kij is the inverse effective mass matrix associated 
with the spin-orbit interaction. The anisotropic effec- 



tive mass can be defined as ^ = Iji^i 



where 



of X, 



■MiijQricEj, 

m + ) ■ 1^ off-diagonal components 

.y vanish in some system, then the diagonal compo- 
nents of the effective mass are equal to the free effective 
mass. This implies that the transverse effective mass is 
zero, and thus, the non-zero transverse SOIF is forbid- 
den in this system. However, for the Rashba-Dresselhaus 
system, we find that the transverse SOIF does not van- 
ish. We will discuss the Rashba-Dresselhaus system in 
the next section. 

We note that the mean SOIF (F^^^) depends on the 
second derivative of the energy dispersion. Unlike the 
force hso] that is related to the conventional 

definition of spin current in k- linear systems, the mean 
SOIF (F'^'^-'^) is generally independent of the definition of 
spin current. We also note that the mean SOIF (F^'~'^) 
is analogous to the semiclassical force acting on a band 
electron (or hole), which can be written as = qE 
and V = 
we have 



dE. 



^ (see Ref. Q). From these two equations. 



dhk 

^ = m^mt = iMfk-lEj; this is in agree- 
ment with the result of the linear response calculation. 
However, it must be pointed out that the force generated 
by spin-orbit interaction has not yet been studied. An 
important indication of our result is that the inclusion 
of the spin-orbit interaction does not change the form of 
the semiclassical force. 
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We now turn to the discussion of the k-hnear system 
where — 0. We note that in the k-hnear system (/iso 

contains only k-hnear terms), we must have qj^^q^. — 

and thus = 0. This means that (F^) has no extra 
term that can cancel (F*^). However, in this case, the 
identity Eq. (fTH) shows that 







9^A„(k) qm 
dkidkj 



(21) 



and thus, Eq. (|17p implies that 



(22) 



In general, 2Re{nk\-^[^^^,hso]\nhy does not equal 
^dt'dk^ ^£j, as can be seen from Eq. ((T7)) . This impor- 
tant relation [Eq. ([22]) ] shows that (F*^) still gives the 
mean SOIF (F'^'^''^) in a k-linear system. Therefore, the 
total force in a k-linear system is still described by Eq. 
([20)1 . The validity of Eq. ([22)l can be further examined 
in the discussion of the Rashba-Dresselhaus system by 
the calculation of Kij [Eq. P^ ]. We will return to this 
examination in the next section. 

In the presence of an electric field, the force experi- 
enced by a band electron (or hole) acts along the di- 
rection of the applied electric field. However, when the 
band structure is anisotropic, the force in the system 
may be perpendicular to the applied electric field even 
in the absence of a magnetic field. The existence of 
a transverse force can be confirmed by calculating the 
off-diagonal component of Kij [Eq. p^ ]. For conve- 
nience, we assume that the external electric field is ap- 
plied in the -|-y direction (or [010] direction), and thus, 
the mean transverse SOIF is {F^'-'^) = qK^yEy and the 
mean longitudinal SOIF is {F^^^) = q{nc -t- Kyy)Ey. 



dkidkj ' 



and 



fc[^(i^)] sin(/)COS(/) and the inte- 



The integrand of Eq. (fTO)) is of the form 
we have the symmetric property: Kij = Kji. If the 
energy dispersion is isotropic, such as in the Rashba 
(A = ak) [2^, Dresselhaus (A = (3k) k-cubic 
Rashba (A = Ofifc^) and wurtzite {aok + (iok^) ^ 
systems, the quantity K^y (and Kyx) always vanishes 
because 

gration of sin(/)COS</> from to 2ti is zero. For the 
spherical Luttinger Hamiltonian = 272 f;^ and 

^±3/2 = "272^^), all off-diagonal components of Kij 
are zero because the dispersion in this system is also 
isotropic. Therefore, the non-vanishing transverse force 
can exist only in an anisotropic system. 

However, it must be pointed out that not all the 
anisotropic systems would have the non-vanishing mean 
transverse SOIF. For the 2D Dresselhaus type system 
along the [110] direction [24], the energy dispersion is 



\/2pkx, and the second derivative of this energy disper- 
sion leads to zero. For the 2D k-cubic Dresselhaus sys- 
tem [1^, [1^, the mean transverse SOIF also vanishes. 
The energy dispersion of the 2D k-cubic Dresselhaus sys- 
tem is A = Pijkkxky, and we find that the band struc- 
ture with n ^ + can be obtained from the other one 
with n = — by using the tt/2 coordinate rotation, i.e., 
kx — ky and ky — > —k^, and vice versa. On the Fermi 

surface, we have ep — ^ 2m'^^ ~ "i^f^okn {4>) sin (f) cos 0, and 
the right-hand side of this equality must be invariant un- 
der (j) ^ —(j>. Using the fact that the function sin cos 

tv'k^ (-4>) 

is an odd function of (/>, we obtain ep = — ^ — + 

)fci(— (?!)) sin(?!)cos(?!) — ^ ^^z}'^^ + /3£)fc-(0) sin0cos0. 



2m 



k^ {(/)). Furthermore, 



Pok'K—ct)) sin(/)cos(/) 
and thus, we have fc+(- 

Eq. p9|) can be rewritten as J^^d4>{k'l — k'L){2 — 
cos^ (/) sin^ (/)) , and thus, it vanishes by using the result 
kn{—4') = fc-n ('/')■ The same conclusion can also be ob- 
tained by means of symmetry consideration. Because the 
energy dispersion is invariant under 7r/2 rotation around 



the z-axis, we obtain K^n^ 



Kyy and 



-Kyx by 



means of K{^) = U{l)KW{f), where t/(f ) is the 2x2 
rotation matrix around z-axis by 7r/2. Due to the result 
that K^y = Kyx, we have K,j.y = 0. This implies that the 
mean transverse SOIF vanishes in the 2D k-cubic Dres- 
selhaus system and charge imbalance does not occur in 
this case. 



The SOIF vector distribution at a 2D fc-plane is shown 
in Fig. 1, where only the transverse SOIF is illustrated 
in the figure. The electric field is applied in the +y 
direction (i.e., [010] direction), and thus, the energy 
dispersion would be displaced toward +y direction in 
hole system and toward —y direction in electron system. 



pSOI 



\Ey\ represents the transverse SOIF 



of the n-th band. The SOIF vector distribution in the 
hole isotropic system is given in (a) and that in the elec- 
tron system is given in (b). The transverse SOIF vector 
of the 2D k-cubic Dresselhaus system is shown in Fig. 1. 
(c). At the angles 0, f , tt, and the SOIF vanishes in 
the isotropic and anisotropic systems. We note that in 
the isotropic system, the electron (hole) concentration in 
the fc-plane is symmetric along the y-axis, and the mean 
transverse SOIF is zero. This implies that the applied 
electric field would not result in the charge imbalance at 
the edges of the sample. In the Rashba-Dresselhaus sys- 
tem ((d) in Fig. 1.), the electron concentrations on the 
left-hand and right-hand sides of the fc-plane are equal. 
However, the mean transverse SOIF does not vanish (see 
the following section). In this system, the non- vanishing 
mean transverse SOIF acts as the driving force and re- 
sults in charge imbalance at the edges of the sample. 
This implies that a non-zero Hall voltage can be detected 
in the Rashba-dresselhaus system when only an electric 
field is present. This is rather different from the classical 
charge Hall effect, where an external magnetic field must 
be applied to the system. 
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FIG. 



1: (Color online) The transverse SOIF vector {F^^°^ = 
'Ey\) distribution. The arrow (red) represents the 



transverse SOIF vector, and the dashed line (blue) represents 
the deviation in the energy band due to the electric field ap- 
plied in the +y direction ([010] direction), (a) hole system 
with isotropic dispersion; (b) electron system with isotropic 
dispersion; (c) 2D k-cubic Dresselhaus system; (d) Rashba 
and Dresselhaus system with a (3; (e) Rashba and Dressel- 
haus system with a = p. 



IV. RASHBA-DRESSELHAUS SYSTEM 



In this section, we study the mean SOIF in the Rashba- 
Dresselhaus system. The Hamiltonian can be written 
^0 = + O'i'^xky - dyk^) - (3{axkx - ayky). 

The energy dispersion is Enk — nA(k), where 

A(k) = ^(q;2 +/32)fc2 -4a/3A:^fcj, = fc7(</>), and 7(0) = 

In our notation, A„(k) = 



V'a^ - 2a/3sin(20). 
— nA and hso = OL{axky - 



(Jykx) 



P{(Jxkx 



A. Non-vanishing longitudinal and transverse SOIF 

In the Rashba-Dresselhaus system, the longitudinal co- 
efficient Kyy can be evaluated analytically: 



■0' 



(23) 



We note that Kyy is always negative in the Rashba- 
Dresselhaus system when a ^ (3. This means that the 
longitudinal SOIF acts in a direction opposite to that of 
Lorentz force. The longitudinal coefficient Kyy is approx- 
imately 10~^(1/A^) in semiconductor with large built-in 
electric field [a - IQ-^eVA) [H, which is very small in 
comparison with the value of the Lorentz force (free car- 
rier density is approximately 10~^(1/A^) in the Rashba- 
Dresselhaus system, see Ref. [11]). The mean longitudi- 
nal SOIF is non-zero in the Rashba system (/? = 0) as 
well as in the Dresselhaus system (a = 0). On the other 
hand, the integrand of the off-diagonal component of Kij 

(transverse coefficient) is ^ = — ^" kxky. The 
transverse coefficient K^^y can be evaluated analytically. 
We find that Kxy also depends on the Rashba and Dres- 
selhaus couplings, and in general, it does not vanish if 



K^ 



47r2?i2 



dSk 



-\a 

a 

a, 2 



d kx 9 ky 



(24) 



Therefore, the system has a net SOIF in the -x di- 
rection, namely, (F^'-'^) = ~\e\KxyEy and Kxy > 
(see Fig.l (d)). When the Rashba coupling is equal 
to the Dresselhaus coupling, the transverse SOIF van- 
ishes, and further the energy dispersion is not isotropic 
(A(a = /3) = V2a\kx - ky\; see Fig. 1. (e)). 

In the Rashba-Dresselhaus system, we have 
d^hso/dkidkj — 0, and the relation in Eq. ([22]) 
follows. This can be examined by the direct calculation 
of the conventional spin Hall conductivity. The force 

is shown to be proportional to the conventional spin 
longitudinal conductivity ayy It can be shown 

(a^ - /32)/^^^d0^Hl^. We have 

for /32 > a\ 



that : 
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^ > P"^ and (T, 

0\ _ qm'^ P 



We obtain 

and (PO) = . 
in agreement with the result of Eq. 



qm 

•2 



yy 
2 



£(a 
a')E, 



Stt P 

P^)Ey for 



a 



for /32 > 



> /32 

a^; this is 
Therefore, 



(Jyky) 



physically the linear response of F*^ to the weak electric 
field is identical to the semiclassical force acting on a 
band electron in this case. Due to the non-vanishing 
transverse force in the Rashba-Dresselhaus system, the 
charge imbalance would exist in the system even in the 
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absence of an external magnetic field because the mean 
transverse SOIF acts as the driving force. 

On the other hand, when an electric field is applied in 
the [110] direction denoted as k'y direction, the effective 



Hamiltonian can be rewritten as 



cry){a — f3)k'^ + -^{(^x + cry){ct + P)k'y, where the energy 

dispersion is "/'{(!>) — \J {ct — cos^ (f) + {a + PY sin^ 0. 
In this case, it can be shown that Kx'y' integral van- 
ishes. It is interesting to note that the conventional 
spin-Hall conductivity in this case is not an universal 
constant [13] ■ We find that the conventional spin-Hall 
conductivity depends on the strength of spin-orbit cou- 
Pling: = + D' « > /5 and al,y, = -^{l + f ), 

a < [i. Furthermore, we find that cr^,^, = — q)j 

Q > /3 and dy,^, = ^(1 — ^), /? > a, from which we 
have an unsymmetrical result cr^iy, ^ ~'^y>x'- This is be- 
cause [110] and [110] are nonequivalent axes in the sense 
that_^the Rashba-Dresselhaus system has the C2v symme- 
try [28]. The resulting spin splittings along the [110] and 
[110] directions are equal to 2{a + f3)k' and 2(a — I3)k' , 
respectively. 



B. The Hall voltage 

When an external in-plane electric field is applied to 
the Rashba-Dresselhaus system, the charge accumula- 
tion generated by the spin Hall current does not result 
in charge imbalance because the system has equal pop- 
ulations of electrons with -\-z and —z polarized spins. 
However, as shown in the above subsection, we find that 
the mean transverse SOIF which is independent of spin 
components does not vanish in the Rashba-Dresselhaus 
system. As a result, the existence of the non-zero trans- 
verse force leads to a charge imbalance between two edges 
of the sample. Furthermore, we can estimate the result- 
ing Hall voltage from the classical point of view. When 
the system reaches an equilibrium, the transverse electric 
field generated by the charge imbalance can balance the 
mean single-particle transverse SOIF which is defined as 



(25) 



where ^ ^ + {^f{o? + /J^)] is the electron 

concentration. Therefore, an equilibrium state requires 
{{F^'^^)) = qEx = q^, where Vh is the Hall voltage, 
and W is the width of a device. On the other hand, 
{{Fi"')) = {K^v/nc)qEy = {K,y/n,)qy^, where Vl is 
the longitudinal voltage and L is the length of a device. 
Therefore, the Hall voltage can be written as 



(26) 



Unlike the classical charge-Hall effect, the non-zero Hall 
voltage occurs in the absence of an external magnetic 



field. Figure 2(a) shows as a function of the Rashba 
coupling, where the Dresselhaus coupling is fixed to (3 = 
7.5 X lQ~'^eVA. As shown in Fig. 2(a), the maximum 
value of Kxy/ric for the Fermi energy being smaller than 
10~^ is about 10"^. We assume that the length-to-width 
ratio is equal to 8 and the longitudinal voltage is about 
10-3 (Volt) [29]. We find that the HaU vohage is nearly 
equal to . 

Figure 1(d) shows the direction of the transverse SOIF 
vector in different k points. If the mean transverse SOIF 
in —X direction is denoted as a minus sign and that in 
+x direction is denoted as plus sign, then the populations 
of electrons corresponding to those different directions of 
the mean SOIF is illustrated in the inset of Fig. 2(b). 
We define a deviation of electron concentration as 



5n = n+ — n_ 



(27) 



where n-|- and n_ refer to the electron concentration with 
the transverse SOIF vector in +x and —x directions, re- 
spectively. In other words, n_ is the sum over all shaded 
regions shown in the inset of Fig. 2(b), and the sum over 
all unshaded regions gives n_|_. Each region ranges from 
zero to the corresponding wave vector, and it follows that 




(28) 



k^<j))d<j)+ / ki{<j))d<j) 



where fc±(0) = ±^7^ + + '-^f 



Equation ((27|) is interesting because it is the origin of 
charge imbalance that requires the population of elec- 
trons at two edges of the sample to be different. In this 
case, we find that the deviation dn/uc grows as the max- 
imum value of Kxy/ric denoted as {Kxy/nc)max increases 
(see Fig. 2(b)). This result explains that the charge im- 
balance can be characterized through the calculation of 
Kxy and clarifies the significance of our expression for the 
mean transverse SOIF. Quantitatively, the calculation of 
Kij [Eq. P^ ] is important because the charge imbalance 
cannot reach to the value 5n, and thus, the resulting Hall 
voltage cannot be determined from 5n. 

We close this subsection by some discussions on gauge 
transformation and the charge Hall current. By using 
Eqs. (|lip and the charge current Ji = qvi in the 

ballistic region can be written as [soj 



hv 



where 



in) 



dA 



in) 



dA, 



(n) 



dk. 



dki 



(29) 



(30) 



is the field strength of n-th band in fc-space, and it sat- 



isfies the antisymmetric property .7-", 



(n) 



(n) 



The 
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FIG. 2: (a) The dimensionless quantity Kxy/ric as a function 
of the Rashba coupling; (b) The deviation of the electron con- 
centration Sn/ric varies with respect to the maximum value 

of Kxy/Uc. 



Berry vector potential in Eq. ([50)1 is defined as A-""* — 
(7ik|(— |nk). We note that the longitudinal charge 
current always vanishes, as shown in Eq. ()30p . where the 

anti-symmetric property of J-'f^'^ guarantees that ^f"^ is 
zero, and this result is independent of any gauge choice. 
Physically, this means that there is no steady longitu- 
dinal current without considering the collision with im- 
purity. In that sense, the existence of a steady charge- 
Hall current generated by the spin Hall effect is possible 
because there is no transverse electric field. However, 
the charge Hall current can not be uniquely determined 
from Eq. (^5]) because the singular gauge transformation 
changes the magnitude of the c harg e Hall conductivity in 
2D spin-orbit coupled systems ^0|. On the other hand, 
the system has equal populations of electrons with -\-z 
and —z polarized spins. Therefore, in Ref. ^SOj], a specific 
gauge is chosen such that the steady charge-Hall current 
is zero. 

Nevertheless, the mean transverse SOIF does not van- 



ishes in the ballistic region, and thus, the non-steady 
charge Hall current exists in the Rashba-Dresselhaus sys- 
tem. As a result, the charge imbalance between two edges 
of the sample can occur in the system. The system can 
still have an equilibrium state in the sense that the trans- 
verse electric field generated by charge imbalance can bal- 
ance the mean single-particle transverse SOIF. 



C. Disorder effect 

In this section, we will consider the weak disorder ef- 
fect. We will show that the transverse SOIF vanishes if 
the calculation is carried out to the order of — in the 

sense that the Fermi energy [ep = h'^kp/2m) is much 
larger than the spin-splitting {Ap = 2'-f{(t))kp). In the 
Rashba-Dresselhaus system, the transverse force is de- 
termined by F*^ because = 0. It can be represented 
by the term proportional to spin current Jy = ^{vy, s^}, 

i.e., = ^{a^ - [3^)J^. Therefore the calculation of 
force is simplified to the calculation of conventional spin- 
Hall conductivity. Because is proportional to J^, we 
have to calculate the spin-longitudinal conductivity tr^^. 

The impurity is assumed to have a (5-potential V(x) = 
VblX]j'^(^ ~ -f^j)' where is the location of the i-ih. 
impurity and Vo is the potential strength. The relax- 
ation time r can be obtained from the calculation of 
the self-energy by using the definition — 2|/toI]|. For 
the impurity potential V^(x) = VblX]i'^(^ ~ the 
self-energy can be calculated by the the following equa- 
tion P: E(e) = i!^^^/dkGfe,«(e), where Gfe,,(e) = 
l/[e — tks — S(e)], and rii is the impurity concentration. 
In the Born approximation [Slj , 



- =^(n,Fo^)(7V+ + 7V_), 

T 



(31) 



where Ns is the density of states of band s. In the 
Rashba-Dresselhaus system, the total density of states 
is {N+ + A^-) = + fc-^) - tI^. and thus, 

{riiVQ) — —. The retarded and advanced Green's func- 
tions are diagonal in the helicity space, that is. 







where 



Gf(e) 



1 



e - Eks + iT 
1 



(32) 

(33a) 
(33b) 



r is defined as T = The spin longitudinal conduc- 
tivity ayy can be evaluated by calculating the following 
Green's function 1321: 



= dkTr[j;G%)VyG^{e% 



(34) 
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where e' = e — Hu! and Vy includes the vertex correction 
[see Eq. ([551) ]. The notation Tr represents the trace in 
the hehcity space s. Here, the electric field is applied in 
the y-direction, and we calculate the spin current in the y 
direction, which has a spin component in the z direction. 
The spin longitudinal conductivity cr*^ can be evaluated 
by using the formula [12] 



cTll= - I de 



e«^(e,e-fiu;), (35) 



where /(e) is the Fermi-Dirac distribution. In the ladder 
approximation, the vertex correction of electric velocity 
Vx satisfies the self-consistent vertex equation js^ : 



'y — + {niV^ ) 



(36) 



The electron velocity Vy can be divided into two parts, 
where Vy is the electron velocity in the 



absence of impurities and Vy is the vertex correction in 
the presence of impurities. The vertex correction of ve- 
locity, Wj,, can be written as Vy = X]ji=o '^m'^'m (^^^ ^^^^ 
Eq.^^7])), and the general solution of Eq. ([55]) is given 



Vy-: 



^c^ ; ^ = 0, 1,2,3, . 



(37) 



The matrix ctq is the two by two identity matrix, and Oi 
(i — 1,2,3 corresponds to x, y, z) are the Pauli matrices. 
We substitute Eq. into Eq. and obtain the 

matrix equation 



Q^LvCv = -Dp ; = 0, 1, 2, 3, 



(38) 



where 



---2Tr[G«a,G-4a^], (39) 



(2^) 

where it can be shown that Q30 ~ Q03 — and 



Df, = (n^Vr 



0') J j^TriC'vyG^a, 



(40) 



It can be shown that Q13, Q23 (and thus Q31 and Q32), 
and 1)3 represent purely interband transitions. If we as- 
sume that the scattering effect does not cause any in- 
terband transition (i.e. this assumption is equivalent to 
Ap >> huj), we can neglect these terms in the calcula- 
tion. Then in this case, the coefficient Cz is zero. On the 
other hand, if we restrict each Qij and Dj to the order 
of Ap/ep^ we can neglect Qqi, Q02 (as well as Qiq and 
Q20), and Do from which we can conclude that cq — 0. 
Therefore, Eq. simplifies to the following two by two 
matrix equation: 



]ii Q12 
I21 Q22 



Di 
D2 



(41) 



In the DC limit w — > (taken before the limit r cx)), 
we obtain Di = f (/^ _ 2) - f /z, D2 = -f ^ + f (/a - 2), 
Qii =2-/1, Q22 = 2-/3, and Q12 = Q21 = h, where 
fi {i — 1, 2, 3) are the following integrals: 



/i 

/2 
f3 



!_ (a sin 0-/3 cos 0)^ 

Jo 

1 /■^'^ ^^(asin</i — /3cos0)(a cost 

Jo 7('^)^ 



/3 sin ( 



1 



27r 



{a cos 4> — P sin ( 



(42) 



here, 7((/') = + /S^ — 2a/3sin(20). After integration, 
we have /i = /a = 1 and /2 = — ^ for a > and 
/2 = — ^ for a < /?. Substituting /i, /2, and into Eq. 
(I4ip . we obtain Cj. and Cj^: 



0) 
0) 



—a/h \ 

~p/h ) 



a ^(3. 



(43) 



However, we have Vy = Vy + CxCx + CyGy = ''^i/ + 

(ca; + %)<^x + (cj^ + i)o'i/j whcrc Vy = hky/m. There- 
fore, the vertex correction cancels the anomalous veloc- 
ity jax + ^o-y Finally, the spin longitudinal conductiv- 
ity (Ty^ is proportional to the term VyVyTrla^G^G^] — 
VyVyTr[azG^aoG^]. It can be shown that the trace di- 
rectly equals zero, that is, Tr[azG^aoG"^] — 0, and thus, 
the spin longitudinal conductivity is zero. This implies 
that the transverse force in the Rashba-Dresselhaus sys- 
tem vanishes when the disorder effect is considered. In 
the same approximation, the spin-Hall conductivity also 
vanishes. This can be seen as follows. The spin-Hall cur- 
rent is proportional to kxCz, and the spin Hall conductiv- 
ity a^y is proportional to the term Vyv';^Tr[azG^aoG^] 
in which the trace term is the same as the spin longitudi- 
nal conductivity, and thus, we obtain the vanishing spin- 
Hall conductivity. The result obtained from Eq. 
agrees with the expected result for the pure Rashba sys - 
tem [32,, j34, j35„ j36] and pure Dresselhaus system [Tsl . 
The influence of impurities on spin-Hall transport in the 
Rashba and Dresselhaus system also has been studied 
in Refs. [IJ and [111 by means of Boltzmann equation 
and Kubo formula, respectively. However, the vertex 
correction [Eq. which is non-zero in the diffusive 

regime [ssl |36| is not taken into account in both calcula- 
tions. 

In short, Eq. (j20p is valid in any spin-orbit coupled 
system with the Hamiltonian given by Eq. The 
charge imbalance would occur in the Rashba-Dresselhaus 
system in the clean limit when the the electric field is 
applied in the system. However, it is shown that the 
disorder can cancel the transverse force, and thus, the 
charge imbalance would not occur in this case. 
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V. CONCLUSIONS 

In conclusion, we show that the mean SOIF (spin-orbit 
interaction induced force, Eq. (|18p ) in any spin-orbit 
coupled systems is determined by the second derivative 
of energy dispersion [Eq. (jl9p ] only. The mean trans- 
verse SOIF acts as the driving force and thus results in 
transverse charge imbalance in the spin-orbit coupled sys- 
tem. It has been shown that the mean transverse SOIF 
vanishes in rotationally invariant systems such as Lut- 
tinger, k-linear Rashba, k-cubic Rashba, k-linear Dressel- 
haus, and wurtzite systems. Furthermore, we find that 
the mean SOIF also vanishes in the 2D k-cubic Dressel- 
haus system, where the energy dispersion is anisotropic. 
Nonetheless, we find that the mean transverse SOIF does 
not vanish in the Rashba-Dresselhaus system. This re- 
sult can be verified by measuring the Hall voltage in the 
absence of an external magnetic field. The estimated 
magnitude of Hall voltage is nearly equal to 1/iV when 
a longitudinal voltage of ImeV is applied. In the pres- 
ence of weak disorder, the SOIF vanishes in the Rashba- 



Dresselhaus system because the anomalous velocity and 
vertex correction accidently cancel each other. 

The electrical measurement of intrinsic spin-Hall effect 
in the ballistic region has recently been carried out for 
HgTe nanostructures [s^. Therefore, we believe that the 
transverse charge imbalance resulting from the non-zero 
transverse SOIF can be detected in the ballistic region 
by measuring the Hall voltage. Hopefully, our inter- 
esting prediction would stimulate measurements of the 
Hall voltage in such spin-orbit coupled systems with an 
anisotropic dispersion as the Rashba-Dresselhaus system 
in the near future. 
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